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Abstract. In this work we study the hemodynamics in a stented artery connected either to a collateral artery 
or to an aneurysmal sac. The blood flow is driven by the pressure drop. Our aim is to characterize the flow-rate and 
the pressure in the contiguous zone to the main artery: using boundary layer theory we construct a homogenized 
first order approximation with respect to e, the size of the stent's wires. This provides an explicit expression of the 
velocity profile through and along the stent. The profile depends only on the input /output pressure data of the 
problem and some homogenized constant quantities: it is explicit. In the collateral artery this gives the fiow-rate. In 
the case of the aneurysm, it shows that : (i) the zero order pressure inside the sac is equal to the averaged pressure 
along the stent in the main artery, (ii) the presence of the stent inverses the rotation of the vortex. Extending 
the tools set up in [S] I27| we prove rigorously that our asymptotic approximation of velocities and pressures is 
first order accurate with respect to e. We derive then new implicit interface conditions that our approximation 
formally satisfies, generalizing our analysis to other possible geometrical configurations. In the last part we provide 
numerical results that illustrate and validate the theoretical approach. 
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1. Introduction. Atherosclerosis and rupture of aneurysm are letlial patliologies of tiie 
cardio-vascular system. A possible therapy consists in introducing a metallic multi-layered stent 
(see fig. 11.11 right). This device slows down the vortices in the aneurysm and doing so favors 
coagulation of the blood inside the sac. This, in turn, avoids possible rupture of the sac. 

In this study we aim to investigate the fluid-dynamics of blood in the presence of a stent. 
We focus on two precise configurations in this context: (i) a stented artery is connected to the 
collateral artery but the aperture of the latter is partially occluded by the presence of the stent 
(see fig. 11.11 left), (ii) a sacular aneurysm is present behind a stented artery (fig. 11.11 middle). 
From the applicative point of view these two situations are of interest since they represent a dual 



Figure 1.1. A sketch of stented arteries: with a collateral artery (left), an aneurysmal sac (middle) and a 
3D example of a real metallic multi-wired stent (right) 

constraint that a stent should optimize somehow: the grid generated by the wires should be coarse 
enough to provide blood to the collateral arteries (for instance iliac arteries in the aorta), at the 
same time the wires should be close enough to have a real effect in terms of velocity reduction in 
the aneurysm. 
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Multi-layer metallic wired stents seem to satisfy both the constraints at the same time. Al- 
though experimentally exhibited [11 [26], these facts needed a better mathematical understanding. 
We give here results in this sense, setting a common framework for both phenomena in the case 
of the Stokes flow. 

Inspired by homogenization techniques applied to the case of rough boundaries [1] 1221 129| we 
construct a first-order multi-scale approximation of the velocity and the pressure. By averaging, we 
get a first order accurate macroscopic description of the fluid flow. Indeed, we compute an explicit 
expression of the velocity through the fictitious interface supporting the stent and separating the 
main artery from the contiguous zone. This formula only depends on the input data of the problem 
and some homogenized constants obtained solving microscopic cell problems. In the case of the 
aneurysmal sac we show rigorously that the zero order pressure in the sac is constant and averaged 
with respect to the pressure in the main artery, which was not known. Then we show that formally 
this leads also to redefine the problem in a new and implicit way in the domain decomposition 
flavor. Actually we obtain a new set of interface conditions along the fictitious interface: while 
for the normal velocity they look similar to those presented in [TOl [21 [9] , the tangential conditions 
arc new to our knowledge. They express a slip velocity in the main artery (as in [IS]), but a 
discontinuous homotctic relationship between horizontal velocities across the interface of the stent 
(see system (|3.4|) '). Our results concern the steady Stokes equations, as in [2], the same interface 
conditions are valid in the case unsteady Navier-Stokes case. 

From the mathematical point of view this paper introduces several novelties. The case of a 
sieve has been widely studied in a different setting in [TT] [T^l [H [H [3] . In these works, the authors 
considered no-slip obstacles set on a surface with various dimensionalities but with a common 
point: the velocity was completely imposed at the inlet/outlet boundaries of the fluid domain. 
Although this could seem a technicality, it influences drastically the limiting regime of the flow. 
Indeed a complete velocity profile is imposed as a Dirichlet condition at the inlet /outlet of the 
domain, so that the total flow-rate through the sieve remains constant whatever e, the size of the 
obstacles: a resistive term appears as a zeroth order limit in the fluid equations. In the context of 
blood flow such a regime seems hard to reach: experiences show that when the wires are to dense 
no transverse flow crosses the stent. This suggests that through the porous interface, blood flow 
should be driven by a pressure drop more that a flxed flow-rate. 

In this direction, Jager and Mikelic considered a pressure driven fluid in [19] . But they studied 
an interface whose thickness was independent on e, which seemed useless for our purpose : the 
diameter of the wires of the stent are dependent on the radius of the artery where the stent should 
be implanted. It appears natural to consider roughness size that varies wrt e in any direction. 
Moreover in this paper we introduce both a tangential and a transverse flow along and trough 
the stent. Indeed, in the limiting regime considered by Jager and Mikelic ^9\, the velocity is 
zero. Here when the collateral artery or a sac arc completely closed by the stent, we still expect 
a Poiseuille profile in the main artery. 

At a more technical level, this work improves the approach developed in [5] [27] in order to 
correct edge oscillations introduced by periodic boundary layers. At the same time, we give an 
appropriate framework to deal with this problem in the case of Stokes equations. Indeed, due to 
the presence of the obstacles, the divergence operator is singular wrt e, this implies degradation 
of convergence results when lifting the non-free divergence terms and estimating the pressure. In 
this frame, we decompose the corrections of the superfluous boundary layer oscillations in two 
parts : 

• on the microscopic side we use weighted Sobolev spaces to describe the behaviour at 
infinity of the vertical corner correctors, defined on a half plane. This provides accurate 
decay rates with respect to e at the macroscopic level near the corner. Indeed, using onto 
mappings between weighted Sobolev spaces we improve decay estimates already derived in 
the scalar case in [51 [27]. Then in the spirit of [2] we construct a microscopic lifting operator 
that allows the vertical correctors to fullfil the Dirichlet condition on the obstacles, 

• a complementary macroscopic corrector is added in a second step, that handles exponen- 
tially decreasing errors far from the corners. 

An attempt to break the periodicity at the inlet/outlet of the domain was done in [20] by using a 
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vertical corrector localized in a tiny strip near the vertical interface. But, decay estimates claimed 
in formula (77) p. 1123 [20] seem to work, to our knowledge, only for a priori estimates of the 
error and are not accurate enough to be used in the very weak estimates. 

We underline as well that in the literature [201 EU [221 [23] error estimates between the direct 
rough solution and the approximations constructed thanks to boundary layer arguments concerned 
the norm of the velocity. In this paper we provide error estimates of the same order for the 
pressure as well in the negative Sobolev norm. This is obtained using the microscopic nature 
of the pressure correctors and in particular thanks to the very precise control of lateral correctors. 
We stress that these vertical correctors play a crucial part in our error analysis at several steps of 
this work. 

The paper is organized as follows: in the two next sections, after some basic notations and 
definitions, we give a detailed review of the results obtained either in the case of a collateral artery 
or a sacular aneurysm. We give in section [4] the abstract results that are used in section [5] in 
order to prove the claims. We provide numerical results showing a first order accuracy also in the 
discrete case in section [6] In Appendix [A] we give proofs of existence, uniqueness and a priori 
estimates for vertical correctors in the weighted Sobolev spaces, while in Appendix [B] we detail 
the results claimed for the periodic boundary layers throughout the paper. 

2. Geometry and problem settings. 

2.1. Geometry. In this study we consider two space dimensions. Let us define by Js one 
or more solid obstacles included in J :=]0, 1[^ of Lipschitz boundaries denoted P in the sense of 
the definition p. 13-14 of Chap. 1 in [28]. We denote by Jf :=]0, \ \^\Js the complementary fluid 
part of iTs in ]0, 1[^. Also, we consider a smooth surface 7a/ strictly contained in J7s and enclosing 
P and we denote the domain contained between 7m and P. Then we define: 
(i) Macroscopic domains: 

The e-pcriodic repetition of J7/ is denoted by and reads: 



The spatial variable giving the position of a point in domains above is a vector called x. 
(ii) The microscopic cell domain: 

As the problem contains a solid interface surrounded by a fluid, the microscopic cell problems 
are set on an infinite strip Z defined as follows 



:= U71oe((i, 0) + J/), where m :== - 



the real e is always chosen such that m is an integer. Then we set 



:=]0,1P, 
n[ :=]0,l[x]e,l[, 
rJi,, := n\ U (]0,l[x{e})U£, 
^2 :HO,l[x]-l,0[, 
Fo :=]0,l[x{0}, 

:=: 17i U To u rj2, 
fie :=f2i,eUroUl72, 

vt' := n[ u rja, 





ri„:={0}x]0,l[, 
Touta :={l}x]0,l[, 
rout,2 :=]0,l[x{-l}, 
Ti :=]0,l[x{l}. 



Z- :=]0,l[x] - cx),0[, 
E:=]0,l[x{0}, 
Z+ :=]0,l[xR+\ J,, 
Z := Z+ U EU Z", 

Zy^^ := Zn]0, l[x]7, i^l (7, e R2 s.t. 7 < I'. 



The microscopic position variable is denoted by y := x/e. 
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(iii) The "corner" microscopic domain: 

In order to handle periodic perturbations on the lateral boundaries rinUr2Urout,i one needs 
to define a microscopic zoom near the corners O := (0,0) and x :— (1,0) of fl^. This leads 
to set the half-plane 11 and the corresponding boundaries as 

n := M+ X M 

N {0}x]0,+oo[, 

D := {0}x] - oo,0[. 

If we choose the obstacle J's to be a single disk, then a graphical illustration depicts the definitions 
above in fig. 12.11 for e = 1/11. 




Tout,! Tin 




fout,! 




Figure 2.1. The macroscopic domains (left) and f! (middle) and the microscopic infinite strip Z (right) 



The exterior normal vector to any domain is denoted by n, if not stated explicitly n is orien- 
tated from Oi towards ^2 on the fictitious interface Fq. The tangent vector is defined as r. 

2.2. Notations and definitions. 

(i) Any two-dimensional vector is denoted by a bold symbol: u := (ui,U2), and single com- 
ponents are scalar and are not bold. The same holds for the function spaces these vectors 
belong to: bold letters denote vector spaces, for instance L^(ri) (i^(ri))^. 

(ii) If 7/ e i?ioc(^) then wc set 

^(y2) / ri{yi,y2)dyi, ?/2 e M, 
Jo 

to be the horizontal average of a function defined on the infinite periodic strip Z . Moreover 
by the double bar we denote a pieccwisc constant function defined on Z as 

f (y) := ?7(+oo)iz+ [y) + '?(-oo)iz- (y), y ^ z, 

whenever the function r/(-) admits finite limits when |j/2| — > oo. We need the values of the 
above function near the origin, thus we set also: 

(iii) For any pair (u,p) € L^(J7) x H^^{VL) we denote by (Tu,p the 2x2 distributional matrix 
reading 

^\i,p{x) ■= Vu — pld2, a.e. a; G fi, 

where Id2 is the identity matrix in M?. The tensor du.p looks like the stress tensor but it is 
not symmetric. This is due to the incompressibility constraint: the Stokes problem can still 
be put in the divergence form with the definition of ciu.p above. 
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(iv) The brackets [•] denote throughout the whole paper the jump of the quantity enclosed across 
fictitious interfaces: across Fq on the macroscopic scale, or across E on the microscopic scale, 
so that for instance 

[ctu,p] := cru,p(a;i,0+) - cru,p(a;i,0"), while :=?/^-^ . 

(v) For every microscopic function rj defined on either Z or H, we denote by 

rie{x) = Tj ^-^ , Vx e fie. 

We also need cut-off functions that wc define here: 

(vi) The cut-off is a scalar function : — > [0, 1] s.t. is a C°°(]R+) monotone decreasing 
function and 

for any positive real z. 

(vii) The "corner" cut-off functions : set ipi := il>{x) and ip2 '■= tpix—x) and ip is a radial monotone 
decreasing cut-off function such that 



Finally set ij{x) := ipiix) + ^2{x). Note that with this definition d^^j = on Fin U Fout,i- 
(viii) The "far from the corner" cut-off function : $ is defined in a complementary manner on 
Tin U Touts U F2 such that 

fV + 3' = l, onFinUFout.iUFa, 
[a„$=0 on FinUFout.i, 

and one shall take for instance $(x) := 1 — ■(/'(O, X2) for all x in fi. 
(ix) For regularity purposes we set Ae to be a cut-off function in the e-neighborhood of the corners 
O and X (p. 1122 [20]). First we set at the microscopic level: 

\{y) := \y2\iB{o,i)iy) + |^i{n\B(o,i)}(y)' '^v ^ n, 

then we define 

Ae(a;) := A (^) M^) + A Mx) + ^x), Vx £ n. 

and an easy computation shows that 

IIAellffi(O) <fc{|log(e)|Ul}, (2.1) 
where the constant k does not depend on e. 
3. Main results. 

3.1. The case of a collateral artery. Wc study the problem : find {u^,p^) solving the 
stationary Stokes equations 

-AUe-|~Vpe=0 in fie, 

div Ue = in fie, 

Pe = Pin on Fin, Pe = PoutS OU Fout,!) Pe = Pout, 2 OU Fout,2, (3.1) 

Ue • T = on Fin U Fouta U Fout,2, 

Ue = onFiUFaUFe, 
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Because of the microscopic structure of , the solution of such a system is complex and expensive 
from the numerical point of view. For this reason throughout this article we use homogenization 
in order to construct approximations of (ue,Pe). This technique decomposes in two steps : 
1) the derivation of a multi-scale asymptotic expansion and the construction of an averaged macro- 
scopic approximation. The first part can be seen as an iterative algorithm with respect to 
powers of e : 

(a) pass to the limit with respect to e and obtain a macroscopic zero order approximation. 
In our case, because of the straight geometry of the main artery and of the boundary 
conditions, the Poiscuillc profile is obtained in f^i and a trivial solution in ^2 : 



Pout,l 




(1 - a;2)a;2eiloi, 



(3.2) 



a:i) +Pout,l2;i)10i +Pout,2la2i 



(b) construct microscopic boundary layers correcting errors made by the zeroth order approx- 
imation on Fj and Fq: we set up in the next section three boundary layers (/3, tt), (T, vd) 
and (XtV) to this purpose. These functions solve periodic microscopic problems (j5.2p . 
((O)) and ((O)) on the strip Z. 

=± =± 

(c) compute the constants that these correctors reach at y2 ~ ±oo: (/3 ,0), (T ,0) and 

= =± =± =± 

[X^T] ). Then subtract them to the correctors. Physically, /3 , T provide a microscopic 
feed-back relative to the horizontal velocity (see the wall-law framework in [521 HD] and 
references therein) whereas the pressure difference represents a microscopic resistivity 
in the flavor of [Hill |6]. 

(d) take into account the homogenized constants on the limit interface Fq by solving a macro- 
scopic problem: find (ui,pi) s.t. 

in Qi U fl2, 
in Qi U fl2, 
on Fi UFz, 




on F 



(3.3) 



N, 



=± 



(xi,0+)/3i 







dx2 . 





ei 



[Po] = 

\m 



on Fo^ 



This macroscopic corrector depends on the zeroth order approximation and the homog- 
enized constants. Due to the explicit form of the Poiseuille profile, the Dirichlet data 
is explicit on both sides of Fq, (nevertheless the solution (ui,pi) is not explicit inside 

U 172). 

(e) go to (jlb[) and correct, on a micrscopic scale, errors made by (ui,pi) on Fe U Fq in order 
to get higher order terms in the asymptotic ansatz. 
2) The second step consists then in averaging this ansatz and obtaining an expansion of the 
macroscopic solutions only. This gives, for instance, at first order : 

U(:(a;) := Uo(a;) -I- eui(x), p^{x) := po{x) + epi{x), Va; G f^i U ri2. 

In particular as • n = e[po]/[^ on Fg, one gets an explicit first order velocity profile across 
Fq. As a consequence, we obtain a new result : 

Proposition 1. The flow-rate in the collateral artery D,2 can be computed explicitly and reads 

Qfo '■= / Uc • ndsi = / [po]dxi ^ -rLr / {pout.l + {Pin ~ Pout.l){l - Xi) ~ Pout,2) dxi 

As stated above [tj] depends only on the geometry of the microscopic obstacle J7s and is 
independent of any other parameter. In the last section of this paper we give some numerical 
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examples that illustrate the accuracy of this result (see fig. 16.51 and l6.6|) . Note that the zeroth 
order approximation does not provide any transverse flow through Fq. Although our results 
provide a first order correction, we underline that in the physiological context the pressures 
(Pin7Pout,i) present in the main artery can be very important compared to Pout.2 : the first 
order flow rate Qpo can thus be quantitatively significant as well. 

In this work we constructed an suitable mathematical framework in order to analyse the error 
made in the two main steps of the construction above. This allows to state the main result of this 
paper: 

Theorem 3.1. There exists a unique pair (ue,pc) G H^(51e) x -^^(^^e) solving problem p.l[) . 
The averaged asymptotic ansatz {ue,p^) belongs to L^(rij) x H~^{flj) for j G {1,2} and satisfies 
the convergence result 



L2(f2iUJ22) 



if-i(f2'iU£,U02) 



< ke 



where | represent any real number strictly less then | and the constant k is independent on e. 

Expressing interface conditions satisfied by (ucp^) on Tq in an implicit way and neglecting 
higher order rests, we show formally that in fact (ue,Pg) solve at first order a new interface 
problem : 



- Aue + Vp, = 
divUf = 
u, = 
u, • r = 

= Pin, on Tin, Pe = Pout, 1 on Tout,!, Pe 



Pout, 2 



dx2 



+ 



u+ • n = u, • n = -^([crn,,pj • n, n) 



in fli U f22, 
in r^i U f22, 
on Ti Ur2, 
on Ttv, 
on Tout, 2, 

on Tn. 



(3.4) 



The horizontal velocity on Fq^ is related to the shear rate trough a kind of mixed boundary 
condition alike to the Beaver, Joseph and Saffeman condition |20] . This implicit relationship 
accounts for the friction effect due to the obstacles that "resist" to the flow in the main artery. 
Nevertheless because the interface separates two domains fli and O2, we obtain a second expression 
between the upper and the lower horizontal velocities ufi and u^i- they are proportional and 
thus discontinuous. To our knowledge this is new. 

On the other hand, the interface condition on the normal velocity could be integrated in the 
Stokes equations as a kind of "strange term" in the spirit of [TOl E] , but as we are at first order with 
respect to e: (i) the strange term is divided by e (in [101 E] this is a zero order term independent 
on e ) (ii) the derivation docs not follow at all the same argumentation. In a forthcoming work we 
study the well-posedness of such a system as well as its consistency with respect to (ucp^) and 
(uc,pe). Because of the particular signs of the homogenized constants but also the discontinuous 
nature of the interface conditions in the tangential direction to Fq, this seems a challenging task. 

3.2. The case of an aneurysm. The framework introduced above can be extended to the 
case of an aneurysm; considering the same domain fig as above we define a new problem : find 
(u£,pe) solving 



- Aue + Vpe = 
div = 

Pe = Pin on Fin, 
Ue • r = 
u, = 



in rig, 
in rig, 

Pe = Pout.i on Fout.l, 
on Fin U Fout.i, 
on FiUF2UFout,2. 



(3.5) 



The main difference resides in the boundary condition imposed on Tout, 2 ■ here we impose a 
complete adherence condition on the velocity ; this closes the output Tout. 2 and transforms the 
collateral artery into an idealized square aneurysm. 

Again, we construct a similar multi-scale asymptotic ansatz. We extract the macroscopic part 
to get a homogenized expansion (ucp^) reading 

Ue{x) := Uo{x) + eui{x), p^{x) -.^ po{x) + epi{x), V.t £ l^i U 172, 

where (uo,po) is again a Poiseuille profile but complemented by an unknown constant pressure Pq 
inside the sac: 

{f \ Pin Poutjl /-. \ Ti 

Po{x)=p+ix)ln,+Poln,, ' ^ " (^.6) 

Pq{x) := pouta + bin - Pout,i)(l - 2:1), *^ 1^ 

Then again (ui,pi) solves a mixed Stokes problem p.3p . the only difference being that Ui = on 
Tout, 2- This gives again a new result: 

Corollary 3.1. The zeroth order pressure is constant in ^2, moreover it satisfies the fol- 
lowing compatibility condition with respect to the pressure in the main artery: 

= PFT / Po{xuO)dxi. 

This gives an explicit velocity profile on Tq which reads: 

u, ■ n ^ -^{p+{xuO) - Po) + Oie^). 



The interface condition exhibited on the normal velocity shows rigorously a phenomenon 
already observed experimentally [HHS]- Set xi^max := niaxa;gro a;i (resp. xi^min := nrin^^gro 2^1) 
and xi := (xi^max + 2:i,min)/2, when xi < xi the pressure jump [po] := Pq{x) — Pq is positive, 
otherwise it is negative. This implies that the first order flow trough the stent is entering fi2 when 
xi < xi and leaving it otherwise. Thus the prosthesis inverses the orientation of the cavitation in 
fl2 with respect to the non-stcnted artery (see fig. 13. ip . 




Figure 3.1. Streamlines and velocity vectors in an aneurysmal sac, with (left) and without a stent (right) 

As stated in the corollary, we show in the next section that in fact the zero order pressure is 
the only constant that insures conservation of mass in fl2- From the medical point of view the 
two claims on pressure and flow are of interest. They quantify and confirm the stabilizing effect 
of a porous stent: besides reducing the stress on the wall of the aneurysm, the stent averages also 
the pressure inside the sac avoiding for instance corner singularities (see fig. 13. 2p . 
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Figure 3.2. Pressure in an aneurysmal sac, with (left) and without a stent (right) 



The geometry presented as an illustration in figures 13.11 and 13.21 does not fit exactly in the 
hypotheses of section [^TH the main difference is the curved circular form of the boundaries of 
Nevertheless the phenomenon observed when is the square ]0,l[x] — 1,0[ still happens when 
^2 has this more physiological shape. 

Again one has a mathematical validation of the formal multi-scale construction 
Theorem 3.2. There exists a unique pair {u^,Pf) <S H^(rie) x L'^{^e) solving (|3.5|) . The 
first order approximation {Ue,p^) belongs to L^(r2i U fl2) x H~^{fli U ^^2)7 moreover we have a 
convergence result that reads 



where | represent any real number strictly less then |, the constant k depends on the data of the 
problem and the domain but not on e. 

We show the same type of result as above : {u^,p^) solve formally the same implicit problem 
()3.4p up to the second order error, but with a homogeneous Dirichlet condition on Tout, 2- 



4. Technical preliminaries. In this section we introduce the basic results that allow to deal 
with the Stokes problem on a perforated domain fij together with the specific boundary conditions 
as in (13.11). 



4.1. Weak solutions for sieve problems. In the spirit of Appendix in [30] we start by 
the definition a restriction operator i?^ acting on functions defined in D, and providing resulting 
functions defined on fi^ and vanishing on T^. 

Definition 4.1. Let V := {v e H^(f2e) s.t. v = onTo and w ■ t ~ Q on T^}, and we 
endow it with the usual norm. The tilde operator refers always to an extension by zero outside 
Of i.e. 



Vu e V,u e H^(r2) s.t. u 



otherwise 



We define the restriction operator G £(H^(r2); H-^(r2e)) s.t. 

(i) u G V implies R^u = u, 

(ii) div u ~ in n implies that div (Rtu) =0 in fie, 

(Hi) There exist three real constants fci,fc2 and independent on e s.t. 



\RM\.2u. ^ < fcl{||u 
Vi?eU||T.2^o > < ^2 I -||U 



L2(o,) :i '^ull"llL2(a) 

|V-ReU||L2(o,) < fc3^||u||Hi(o)- 



+ 1|VU|1^2(0)4 f , 



Lemma 4.2. There exists an opertor Rf^ in the sense of Definition\j7Jj 

Proof. The restriction operator i?^ is constructed exactly as in Lemma 3 and 4 in the Appendix 
by L. Tartar in [30], namely for a given u e H^(J') there exists a unique pair {v,q) S H^^Jm) x 
(i^(JM)/M) satisfying: 



( - Av + Vg = -Au 
div V ~ div u 



\Jm\ Jj. 



V = 
L V = u 



in Jm, 
div u dy, in Ja/ 

on P, 
on 7a/, 



There exists a constant k independent on u s.t. !lv||jji^j-^^j < fc||u|| jjij-j-j. By setting 

Ru{y) { v{y) if y G Jm 
if y e J. 

we evidently have ||-Ru||jjij.jrj < fc||u|| jji^jj . and Ru coincides with u if u = on J'^ and div u = 
implies div (Ru) = 0. 

Now let u e H^(f2). For any given x G fe, we set yi := Xi/e — E{xi/e), y2 '■= X2/€ and 
i := E(xi/e) where E(-) is the lower integer part of its real argument. For each i we define a 
function : — > R s.t. u^{y) := u(a;). This allows us to set R^u as 



u{x) 



if x e n', 



RM^)--={ ^ 

2_^Ru^{xi/e~ i,X2/e)l^(^j+,ei){x) li x e £e 



This definition implies obviously that |li?eu||jji(j2,> = ||u|| jjwq,n and we focus on S^. 



i=0 



The key point of the proof are now estimates on the gradient. Taking a regular function u e 'D{n), 
one obtains in a similar way as above: 



1, 



u(a;) = u(xi, 0) + / u(a;i, s)(is 



Now one writes 



which, after taking the square, integrating on £^ and using Cauchy- Schwartz gives 
thanks to the continuity of the trace operator ||7(u)||^^ ^ < ^i|u||jjijf^j, one has: 



'{6||u||^,(,,)+62||Vu||i.(£^).}, 
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(4.1) 



where the constant k' docs not depend on e. Using this last inequality in (|4.1[) ends the proof. □ 
Definition 4.1. We define the corresponding lifting operator S^u (R^ — Id2)u. For every 
X in £j there exist a unique i := E{x/e) € {0, . . . , -}, yi = Xi/e — i and y2 ~ ^2!^ s.t. ifvi solves 



( - Av, + Vg = 
div V,- = 



in Jm + iei, 
div Uj dy, in Jm + iei 

on P + iei, 
on 7a/ + iei, 



then one sets Seu{x) :— J2i=o^iiy)'^JM+iei for every x G il^- One has estimates similar to those 
of the restriction operator 

ll'5'eu||L2(o,) = \\Seu\\-i^2^c,) < /c|||uj|^2(£^) +e||Vu||^2(^^)|, 

l|V5^u|lL^(f2,) = I|V5,u!|l2(£^) < |^||U||^.(£^) + i|Vuj|^.(£^)| < -^JHv ■ 

Proposition 4.3. Let g e L^i^le) there exists at least one vector v e V s.t. 

k 

divv = g in a, |v|Hi(f2^) < -A9\\L-^(n,y 



Proof. We extend g by zero in \ which we denote g, wc use Lemma III. 3.1 and Theorem 
III. 3.1 in [IS] stating that there exists w £ H^{^) s.t. 

divw = g in O, |v|jji(j^) < A:||.g||i2(f2), 

where the constant k does not depend on e. Using the restriction operator defined in the proof 
of Lemma FOl one sets then 

thanks to the estimates that the restriction operator satisfies, one gets the desired result: 

k k' k' 

□ 

Thanks to the latter proposition one easily gets by duality arguments as in p. 374 in the 
Appendix in [30] 

Proposition 4.4. There exists a constant k independent on e s.t. for every distribution 
p e V'iQe) s.t. Vp e V', one has 



At this stage we can derive existence and uniqueness as well as a priori estimates for the 
solutions of the problem: given {i,g,h) e y' x ^^(rig) x i7~2(r7v) find {u,p) s.t. 



f - Au + Vp = f in fie, 
div u = g in 
u = on Td, 

ur = Ol 
I P = h J 
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(4.2) 



on F 



Theorem 4.5. If the data of problem (|4.2p are s.t. {f,g,h) e V' x i^(rJe) x H^^{Tn) then 
there exists a unique solution (u,p) G V x i^(ile), moreover one has: 



I|u|Ihi(o,) + \\p\\mn') + V^Mlhc,) < |l|f|l v + -^hh^n.) + \\9 - /^ll^- i (p^^ 

where the constant k is independent on e. 

Proof. The existence and uniqueness of (u,p) E V x L^{^1^) arc standard results of the 
literature (sec for instance [13] and references therein) . We focus here on the control of the norms 
for this solution pair. Lifting the divergence source term provides easily a priori estimates on u: 

W^^Wmn.) < {llfllv + ^II.9!Il^(o.) + 11.9 - ^!l^-i(r„)} ' 



Then we split fie and restate problem (j4.2p on il' , having for the pressure that 

-Au + Vp = f in n', 

which gives 



l|Vp|lH-i(jy) < k\\i + ^^\\H-Hn') < k{m\v' + I|u|Ihi(o.) + M 
but on this domain there exists a constant independent on e s.t. 

\\p\\mn') < fc||Vp|lH-i(o')' 

which gives the same error estimates for the gradient of the velocity as well as for the pressure in 
r2'. Unfortunately because of the presence of the obstacles, in the rough layer one has only that 

Ml^c.) < Mlhu,) < ^llVpllv = ^||f + Au||^„ 

which, by using again the a priori estimates of u in fi^, gives the final estimate. □ 

4.2. Very weak solutions. We recall here the framework of "very weak" solutions originally 
introduced in US [H [H]. 

Definition 4.6. Let uj be an open bounded connected domain whose boundary du is split in 
two disjoint parts doJo CLnd Oujn- It is said to satisfy the regularity property x for the 
Stokes problem if for all F G L^(w) and every G £ Hq^lo) the solutions of the problem 

in w, 
in UJ, 

on dujD, (4.3) 
on dujNi 

satisfy T G H^(u;), X G H^iui) and if there exists C'l = Ci(cj) s.t. 

WrWuH.) + PIIhm^) ^ C'l {||f|Il.(.) + IIvg||l.(^)} . 

Following exactly the same proof as in Appendix A in |12j one shows 

Theorem 4.7. If uj satisfies the regularity property of definition \4.6\ above, then there exists 
a unique solution (u,p) G L^(w) x H~^{uj) solving: 

in UJ, 
in UJ, 

on dujN, 

on dujj), 





(4.4) 



^, u = m. 
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provided that the data satisfy: f € V^uj), g E L'^{u!), h E H ^{Oun), £ E L^((9cjjv); G 
H^^{dujpf) and m G L^(9w£)). Moreover there exists C2 = C2{<^) s.t. 

II^IIl^M + M\h-h^) < ^2 {lim|iL2(a^^) + ||f|lv'(<.) + IIsIIl^m 

+ IKIlL2(at^„) \ ■ 

where := {v S H^(cj) s.<. v = on and v • t = on 9wAr} and ^'(w) is its dua/. WTe 

denote by "very weak" solution such a pair (u,p). 

Theorem 4.8. // the pair (u,p) is a weak solution of problem (|4.2p . one /las t/ien t/ie very 
weak estimates: 

WAi^^iV) + \\p\\H-^{n') < {llfllv + \/e||u|lHi(n.) + Il5 - /^llH-i(riv) + II5|Il2(o.) } ' 
where the constant k is independent on e. Moreover in the rough layer one has: 

I1u|Il2(£.) + lbllff-i(£,) < {l|Vu||i2(o,)4 + lbllL2(£^)} , 

where again the generic constant k is independent on e. 

Proof. The first estimate follows by applying Theorem 14.71 with w := O'. It is easy to show 
that actually in this doamin, the constants present in the very weak estimates are independent on 
e: the obstacles are not part of fl' . Thus one has: 

I|u|Il2(o') + WpWh-HW) < '={l|u!lL^({:r2=o}u{x2=e}) + l|fllv'(n') + hh^n') + 11.9 - /^ll_ff-i(r„nafy)} > 
< k' {v^!|Vuj|L2(£^) + i|filv'(oo + hWmn,} + 11.9 - /^llH-i(r«)} ' 

where we used Poincare estimates knowing that u vanishes on T^. It remains to consider the rough 
layer C^. There, we have 

I|u|Il2(£.) + \\p\\h-^c,) < e {||u|Ihi(£.) + bllL^(£,)} ' (4-5) 

where the tl^{Ce) regularity is obtained using Theorem 14. 5 1 Indeed the estimates on the velocity 
come using Poincare estimates at the microscopic level as in Lemma 3.2 in jllj . the pressure 
estimate is obtained by duality: by definition of the dual norm one has 

WpWh-uc,)"^ sup <P,'P>H-Km^ 

where Hl{Ce) denotes the set of functions in H^{Ce) vanishing on FqUJO, l[x{e} U dCe- As p 
belongs to L^{Cf) the duality bracket can be transformed into an integral, namely 

<P,^>H-KHt = P^dx< \\p\\l2[c,)\\v\\l2(c,) < 4p\\L2{cJMm{c,v 

taking the sup over all functions in HQ{Ce), one concludes the norm correspondence. □ 
5. Proof of the main results. 

5.1. The case of a collateral artery. In what follows we set both Pout.i and Pout,2 to be 
zero for simplicity. The results remain valid for any fixed constants Pout.i and Pout,2 as well. 
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di , 



5.1.1. The zero order term. When e goes to zero we show m a first step that {u^,pe) 
converges to (uo,po) the Poiseuille profile stated in p.2|) . which solves in ft^ : 



- Auo + Vpo = [o-uo,po] • n(5ro 
div Uo = 
uo = 

Uo • T = 

Po = Pin on Tin, po = 



in r2e, 
in ilf, 
onTiUTa, 
on Tat, 
on Tout,! U Tout, 2, 

on r,. 



Theorem 5.1. For any fixed e, there exists a unique solution {u^,pe) G V x L^{il^) of the 
problem (|3.ip . Moreover, one has 



where the constant k does not depend on e. On the other hand one can prove that: 

- uo|lL2(n,) + Ibe - PoWh-^q'uc,) ^ 



Proof. Existence and uniqueness of the solutions of problem p.ip conic from the standard 
theory of mixed problems [151 IE] j Theorem 14.51 gives more precisely 



where the constant k is independent on e. As Uq docs not satisfy homogeneous boundary conditions 
we use the restriction operator already presented in the section above. Namely we set: 



these variables solve: 

- Au + Vp = 
div u 
u = 

U T = Ol 

p = 0\ 



u:=Ue-i?eUo, p:^Pe-po, 



A(uo - i?eUo) + [o-uo.po] ■ n(5ro = A(S'eUo) + [ctuo.pq] • n^Po in ft 



in fig, 
on Td, 

on Tff, 



where the lifting operator is given in Definition 14.11 Thanks to Theorem 14. 5[ one has then 
directly: 

I|V"IIl2(J2,) + IIpIIl2(o') + V^WPllmc) ^ ll^-^^uo + [ctuo.po] • n<5rollv' 
Thanks to the vicinity of Fg , one deduces easily some trace inequalities |11| : 

mmro) < V4^nLHn.)2, V* e H\n,) s.t. = on T,, (5.1) 
this estimate allows us to conclude that 

sup / ([ctuo.po] •n,*)dxi < ||[cruo,po] ' "11 L2(ro) H * Hl^ (Fq) ^ \/e||[c^uo,po] • n|lL2(ro)ll*llv 



The specific form of the lifting S^Uq allows to write: 



|A(5,uo) 



|V(5,uo)||^.(^^ 



< 



1, 



Uo r 2 



iVuol 



<kV~e, 
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where we used the exphcit form of the Poiseuille profile in the rough layer. Using Theorem 14.81 
one has then 



< ke. 



One has then easily also that 

||Ue - Uo||l2(£^) < ||u, - i?,Uo| 



\ReUo 



uoilL2(£.) < 



I'5'"oI1l2(£^) < ke. 



Estimates above show a threefold error: the Dirichlet error on T^, the jump of the gradient 
of the velocity in the horizontal direction across To, and the pressure jump across Tq. In order to 
correct these errors we solve three microscopic boundary layer problems. 

5.1.2. The Dirichlet correction. The first boundary layer corrects the Dirichlet error on 
r,;. It is very alike to the one introduced in the wall-laws setting [501 UE]. Namely we solve the 
problem: find (/3,7r) such that 



- A/3 + Vtt 
div /3 = 

/3 = -2/261 

132^0 

{j3, tt) arc 1 - 







in Z, 

in Z, 
on P, 
I2/2I 00, 



(5.2) 



periodic in the yi direction. 



We define as in [TS| p. 56, the homogeneous Sobolev space := {v e 'D'{Z), s.t. Vv G 

(L^(Z))^}. Moreover we denote by D\''^{Z) the subset of functions belonging to D^''^{Z) and 
vanishing on P. 

Proposition 2. There exists a unique solution (/3, tt) g D^''^{Z) x Lf^^{Z), n being defined 
up to a constant. Moreover, one has: 



the convergence being exponential with rate 7^ and 

' ^2(2/2) = 0, Vy2 GM\]0,y2,p[, 

< ^1(2/2) = -\Js\- |V/3|2,(^j, +^i(0), Vy2 > 2/2,p, 

.^1(2/2) =^i(0), Vy2<0, 



±00 



where j/2,p '■= max.ygp 2/2 md |j7s| is the 2d-volume of the obstacle Js For sake of conciseness the 
proof is given in the Appendix [B] 

5.1.3. Shear rate jump correction. The second boundary layer corrects the jump of the 
normal derivative of the axial velocity: we introduce a source term that accounts for a unit jump 
in the horizontal component but on the microscopic scale. Namely, we look for (T,-cc7) solving: 



{ - AT + Vtx7 = feei in Z, 

div T = in Z, 

T = on P, 

T2 ^ |2/2|-5^oo, 
(T,tu) are 1 — periodic in the y\ direction. 



(5.3) 



Again we give some basic results and the behaviour at infinity of this corrector. 

Proposition 3. There exists a unique (Y,vu) G D\''^{Z) x L'^^^{Z), w being defined up to a 
constant. Moreover, one has: 



riy) 



=± 

ei, 
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2/2 ±00, 



T2(2/2)=0 VyaeK, 
Ti(2/2) -Ti(0)+^i(0) Vj/2>y2,p, 
Ti(y2) = Ti(0) = ||VT||^.(^), Vy2<0, 

where 2/2, p := maxygpj/2- 

The reader finds again the proof in Appendix |B] 

5.1.4. The pressure jump. In order to cancel the pressure jump [po], we use a corrector 
similar to the one introduced and widely studied for a flat sieve in p. 25: 

' - Ax + V?7 = in Z, 

div X = in Z, 

X = Q on P, (5.4) 

X2 ^ -1, |?/2|^00 

(x, 77) are 1 — periodic in the j/i direction. 

As in the proof of Proposition [21 one repeats the arguments of Appendix |B] in order to obtain 
similarly to [11] : 

Proposition 4. There exists a unique solution {XiV) S D^''^{Z) x Lf^^{Z) of system (|5.4p . 
77 feemg defined up to a constant. Moreover, one has 



-62, 



the convergence being exponential with rate 7,^ and there exists two constants r]{+oo) and r]{—oo) 
depending only on the geometry of P such that 



One then proves: 



|Vxl 



This corrector will be used in the sequel, but we already utilize it to give a first result on the 
average of tt and tu 

Corollary 5.1. The solutions (/3,7r) and {Y^w) solving respectively (j5.2p and (|5.3[) satisfy : 
7r(y2) = and w{y2) = 0, V?/2 e M-U]?;2.p, +00 [ 

For the proof see again Appendix [BJ 

As explained in Remark [T] below, we need in section lS. 1.61 a higher order corrector that solves 
the problem: find (>f, /i) s.t. 

- A^^ + V^i = -2(Vx - (?? - ^)Id2).ei in Z, 

div >f = in Z, 

>f = on P, 

><2 ^ 1 12/2! ^ 00, 

(>f, /i) are 1 — periodic in the yi direction. 

Proposition 5. There exists a unique solution £ D^'^{Z) x Ly^^{Z), ji being defined 

up to a constant. One has also exponential convergence towards constants with rate 7>c ; 

>f — >■ >?, /i — > /I, when IJ/2I — >■ 00. 
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Moreover one has the relationships between values at 1/2 = ±00 

[P] = -2 J^Xiiv~^))dy, [^1] = -2^(cr^_(,,„=).ei,/3 + y2ei)dy. 

The proof is exactly the same as for Propositions [2] and [3] and thus is left to the reader 
In what follows we use the e-scaling of all boundary layers above, namely we set: 

(3,{x):^f3[^), T,(x-):==t(^), x.(^):=x(f), >c,{x) := >c , 

the same notation holds for pressure terms as well 



5.1.5. Vertical correctors on Fin U Tout,! U Above boundary layers are periodic; their 
oscillations perturb homogeneous Dirichlet as well as Neumann stress boundary conditions on 
Tin U Tout,! U The perturbation on these boundaries is 0(1), due to the vicinity of these edges 
to the geometrical perturbation T^. We introduce vertical boundary correctors defined on a half- 
plane n. Each of them accounts for perturbations induced by the periodic boundary layers on 
Fin U Tout,! U r2 in the very vicinity of corners O and x. These correctors solve at the microscopic 
scale the problems: 



'/3 



in n, ( AwT + S70r =0 in H, r - Aw^ + = in H, 

in n, div = in E, 

on £>, S = -(x - f A) on D, 
on N, = -(X2 - X2A) on N, 

on N, [dx = -iv-^) on TV, 

(5.5) 

and {w^, 9^) solves a similar system lifting {xX — x, fi — fi) on DU N. Note that the domain 11 
does not contain any obstacles, we should use a restriction operator on the velocity vectors in 
order to handle this feature (see below ()5.7|) ). We define the usual weighted Sobolev space [TBI [3]. 
for all (to,_p, a) e N X [1, oo[xM: 



- Aw/3 + 
divw^ = 

w^ = -(/3-^A) 

0p = -TT 



in n, 
in n, 

on D, 
on N, 
on N, 



— Awx + V^T = 
div wx = 

WT = -(T -TA) 

WT,2 = — T2 

6r = —fi 



'(n) := e 2?'(n) s.t. G LP(n), < |A| < 771 } 



where p := (1 + It/P)^. We endow this space with the corresponding weighted norm. By density 
arguments one proves that dual spaces of W^'P{IV) are distributions and we set in the rest of this 
work 



W~ 



'(n) := (W^™^f(n))', y{m,p,a) e N X [l,oo[xl 



Here we extend results obtained for mixed boundary conditions and the rough Laplace equation 
in [5J to the case of the Stokes equations. In the appendix we give the extensive proof of the 
crucial claim: 

Theorem 5.2. Thanks to the exponential decrease to zero of the boundary data in (|5.5p . there 
exists a unique solution (wi,0i) G W;|^'^(n)^ x PT/|^'^(n) for i G {/3, T,x, ^^l, for every real a s.t. 
\a\<l. 

Remark 5.1. The weight exponent a provided by this result on the microscopic scale is 
important. It accounts for the behaviour when p goes infinity of the vertical correctors above. The 
decay properties so described are used in Lemma I5.il in order to quantify, in terms of powers 
of e, the impact of the perturbation induced by the periodic correctors on the macroscopic lateral 
Dirichlet and Neumann boundary conditions: the greater a the smaller the error in terms of powers 
of e. So we assume a very close to 1. 



The Poiseuille profile admits an explicit form p.2p and thus its derivative wrt X2 reads 
9x2^0,1 = (Pin — Pout,i)(l ~ '2,X2)/21qi. For the rest of the paper we implicitly assume 9x2^0,1 to 
be evaluated at X2 =0+: it is constant and reads 



dx2 



^"0,1 , „ + . Pin - Pout.l Pin 

^(^^1,0+) = 2^ = ^ 
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(5.6) 



We set for i G {/3, T, k}, 



We,i(a;) := CiR^Wi ( - ) ipi{xi) + CiR^w 



j(a;) := y-j ipixi) + c^O^ ( — - — j ^2(2:), 



(5.7) 



where we used the restriction operator R^ of definition 14.11 while (w,;,^^) solve similar problems 
as (|5.5p but on the halfspace n~ := R_ x R, and the constants c,; (resp. c,) denote 



c/3 := 



dx2 
duo 



-(O), cr 



dx2 
duo,i 
dx2 



(O), c^:=igl(0), 



Pin, 



bo] 

m 



For the particular explicit zero order solution (uo,po) expressed in p.2p . is the only constant 
for which Ci ^ q. As the analysis carried below on Fin U F2 is exactly the same on Fout,i U F2 we 
implicitly assume that when terms appear containing constants c.j, ipi, and 6i similar expressions 
with Ci,ip2,'^i and 9i are considered as well. 

Lemma 5.1. Defining vertical correctors {'We,i,0e,i) with i G {f3,T ,x, one has the esti- 
mates: 

• In the whole domain: 



|A(eWe^.O - V6'e,i||v' < ke^ , Hedivw, 



where the constant k is independent on e, and 1 is any constant strictly less than 1. 
• In U,' one has: 



|A(ew,,,) - V0,,,||H-i(n') < ^e'^", Hedivw, 



< fee 



l+a 



where a is any positive real smaller than 1, and k is again independent on e. 
• On F^v one has 



1 ediv 1 



< fcexp ( — 



• /n f2i U fl2 one has: 

||w,,,||l.(j,^) < fee", J €{1,2}. 
Proof. All easy calculation shows that for any i G {/3, T, >f} one has: 



J := -A(ew,^0 + Ve,^, = -A^ ew, - U' + 6*, - V 



Sf-w,, ( - 



9, (^) Idz) V,V^ - (ew, (^) A,V) + A, (eV^.w, (^)) 



One has then that 



\J\\v'<k 



+ 



A. 



v 



Id2 j V^i/' - i^ewi ) A:r?A 
=: h + h 



L2(n,) 



We split Ji in two parts /i.i and /i.2. The first part is estimates as: 

III < / (|Vw,|2 + < ke^ [ (|Vw,|2 + 0^)^^ 

Jn Jnnli.#-rxio,7r[ 



<e2|||Vw,|+a,||?^o,. 



iVc, (n) 



sup < fcg2(l+a)^ 
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Similarly the second part is estimated as well by 



< ke 



M^rdrdO 



/nn]i,^[x]o,7r[ V P 
Passing then to I2 one has easily that 



p^°'rdrd0 



sup p 



2-2q 



eA ((5,w,) (-) 



v 



((5,w,) (-) ^) 



(5.8) 



Now thanks to the estimates on the lift 
2 



<e1|5,w. 



»llL2(B(0,i)) 



L2(B(0,i))* 



and the same way one gets: 
2 



< k 



Putting together last two estimates in ()5.8p one obtains the first result of the claim. The result 
on the divergence follows the same lines. 

On fl' the result is more straightforward since R^vfj = for all i e {(3, T, x^ Then using 
again the correspondance between macroscopic powers of e and microscopic weighted spaces one 
gets easily the result. 

On Fat one has that 



edivw<:,i = ediv (^w^ V'j = eVV' • (^-j = edr4'i^i ■ t) 

because dnip = on this boundary. On the microscopic scale the support of Ot-iP is located in ^ 
and thus one has 



I ediv 1 



'H 2(^J^ 



< lledivw, 



< e exp 



Then, we define the complete vertical corrector as 
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Vx e fie 



in 

in rif 



on r 



N, 



where (W, S) solve the system of equations on the macroscopic domain fi^ : 

AW + VS-^O, 
divW = 

W • r = 6 {c^(/3, - ^) + cr (T, - T) + c^(x, - f ), +ec^(><. - ^)} • r $ 

S = {cpTT, + c-rrue + c^iVc - ^) + ec^(/^e - P)} * 

W = e [cf3{(3^ -^)+ cr (T, - Y) + c^(x, - ^) + ec^(><. - ^)} $ on Td 

(5.9) 

where one notes that We = on because of the support of <f>. 

Proposition 6. There exists a unique solution (W, S*) G ^.^{ilc) x of system (|5.9p . 

moreover one has: 



IWI 



< fee 



where the exponential rate 7 and i/ie constant k do not depend on e. 
Proof. Setting 



{ dx2 



TZ:= [e 



dx2 



(Y,-T) + ^(X,-X) 



and W := W — 7?.; by the standard theory for mixed problems [El [13], there exists a unique 
solution (W, S) of the lifted problem. One has also a priori estimates : 



W 



Hi(r!.) 



\\S\\mn,_) < llA7^||H-l(^^, 



||div7^| 



l^ll 



(Fj, 



Thanks to the crucial presence of the cut-off function $ and the exponential decrease of rate 
7 := min(7^, 7-r , 7^1 l>c) of all the microscopic correctors, one gets the exponential decrease of the 
rhs in the previous estimates. Because it is also trivial to show that ||7^||Hi(r2 ) — one ends 

the proof. □ 



5.1.6. The complete first order approximation. Having introduced every single element, 
we built a complete first order approximation. We define the full boundary layer corrector: 



Uo 



dx2 

2 J -Pin 



(/3. - /3) + 



9mo,i 



dx2 



(T,-T) + M(X, 



Ul 



-=-(>tfe - >C) +U2 

m 



Po 



I 0X2 



dx2 



bo] 



(5.10) 



roe + -r^iVc - ?7) + epi > + e-s-(Me - P) + f- P2 + 



where the normal derivative 3x2^^.1 is defined in (j5.6p and where the first order and second order 
macroscopic correctors (ui,pi) and (u2,p2) solve respectively p.3p and 

( -Au2 + Vp2 = in f2i U f^2, 
divu2 = in fii U r22. 



U2 =0 

U2-r = 0l 

P2 = J 

Pin=± 

m 



on Td, 
on Tat, 



(5.11) 



on r. 
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Problems (|3.3p and (|5.1ip are defined on two separate domains fti and il.2' two distinct values are 
given as Dirichlet boundary conditions to the horizontal component of the velocity on Fq. This 
is due to the different values of the constants whom the boundary layer correctors f3 and T tend 
to at + and - infinity. Thus the velocity vectors Ui and U2 are not only discontinuous across 
To but also multi-valued at the corners O and x. It appears then clearly that Ui and U2 cannot 
belong to Hi(f]i U 1^2). For this reason, we use the concept of very weak solution introduced in 
the section above. Because ft[ and ^2 are convex polygons in R^, they fulfill regularity conditions 
of definition 14.61 (see example 2.1 p 53 in [l2]). This allows to use Theorem 14.71 in order to obtain 

Corollary 5.2. The pairs of functions (ui,pi) and (u2,p2) solving problems (|3.3p and 
(|5.11|) exist and are unique "very weak" solutions in L^(r2i U ^2) x _ff"^(f2i U il2) ■ 

For some technical reasons appearing later on, one needs to set up an intermediate pair of 
functions (u^jPj) solving : 



- Au^ + Vp^ 
div = 

diiQ. 



0, 



Ui = - 



dx2 



■/3- 



dx2 



[Po] = 



\ 0x2 



dx2 



bo] 



in ill U 0,2 
in VLi U Q.2 

on Fi U F2 



on Fin U Fout.l U Fout,2 



(5.12) 







on Fo 



We define implicitly the same second order problem whose solutions we denote in the same fashion 
(u2 ,^2)- These are regularized versions of problem (|3.3I) and (|5.1ip where we lifted the constants 
from the interface Fq. Indeed the function Ui is multi- valued in the corners and x multiplying 
the specific cut-off function Ac on these corners insures that: 

Proposition 5.3. There exists a unique solution (u^jp^) solving (|5.12[) . Moreover one has 
the estimates: 



ll"i llHi(n) 

Taking the restricition to il^ of one has 



P^\\L^n) - 



where the generic constant k is independent on e. 
Proof. We denote 



lL2(n,)' 



dx2 



m 



On each subdomain by lifting the Dirichlet data one obtains in a classical way 



1 llHi(niua2) 



1^1 llL2(niun2) - ^ 



= k{\\Vg\\L2(mun^2) 
<fc{l|Vg|L2(n,un.) 

<fc{l|Ae|lHl(o)+2|^2ll|l|lL^(0,e)} 



{llg|lH-i(aiun2) + Hivg 
+ l|divg||i2(o,ufi2) 
+ Ildivg|li2(s^^un2) 



L2(aiun2) 



lldiv 



^"H-3(ri„uro„t,i)} 



div ell 1 , } 

divg|lL-(r„ur„„,.i)} 



where the constant k does not depend on e. Note that the latter equality is true since d^Xe 
vanishes on the boundaries of Q and X • ei =0. Now using the estimate on the gradient of 
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Ae (|2.ip . one recovers the first a priori estiraate. Working on ili^e, when writing the system that 
(RciLijPi) solve, one gets easily that 

and because = on To one uses the Poincare inequality on the first term in the rhs above in 
order to obtain; 



On il2, there are no obstacles i.e. R^Ui = Uj which ends the proof. □ 

Remark: 1. In the error estimates developed in the next sections, one applies the momentum 
operator to the term: 

^«Xe-f),^-^)- (5.13) 

Because [pq] depends on x, the rest is not zero: among others a 0(1) double product of gradients 
remains, it reads 

2Vyix. - f ) • ^ = 2^V,(x. - f ) • e, = -2^V,(x. - f ) • er. 
This term could be estimated directly in the L^-norm, giving 

which is a zeroth order error. Needing better error estimates, we add the second order term in the 
asymptotic ansatz (|5.10p that reads: 

~dxiPa{e^{>Ce - e(/x, - p)) 

The Stokes operator applied to this corrector cancels exactly the double product above. The diver- 
gence of the velocity part of (|5.13p gi 
norm, so we should not correct it. 



gence of the velocity part of (|5.13p gives as well a cross term : this is already of order e 2 in the 



5.1.7. A priori estimates. We consider here a complete boundary layer approximation 
containing a regularized macroscopic correctors (u^,p^) and (£12,^2) reading: 



Po 



/9uo,i ^ 



dx2 



dx2 
duo.i 



dx2 



m 
bo] , 



1=] " 
m 

Pin 



(5.14) 



Note that the problems ()5.5p that the vertical correctors We,i solve, account for perturbations 
induced by the periodic boundary layers (3^, Y^, x^, and u^, U2 on Fin U Tout. 1 Ur2, this explains 



the presence of A in the definition of boundary terms in (|5.5p . We recall that these correctors are 
included in the global corrector W^. 

Theorem 5.4. The full boundary layer approximation {U^ ,V^) defined in (|5.14[1 is a first 
order approximation of the exact solution {\i^,p^) i.e. 



U 



Ihi(o, 



+ m 



\fA\Pe 



where the constant k is independent on e. 
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Proof. Set ii -.^ — and p := — , they satisfy: 



- Au + Vj5 = ^ eAwe,j - V8e 

JG{/3,r,x} 

divu= ediv i + O ^£2^ 

i:e{/3.T,x} 



in rif 



u = 







u = 


-( 




u • 


T = 




p = 


= 


1 



9X9 



2:261 + + e^U2 



on TiUTa, 
on Tf , 



on r 



There are two kind of errors: the first is due to the locahsation of vertical correctors and is treated 
thanks to Lemma 15.11 the second is due to the macroscopic approximations that do not satisfy 
the Dirichlet condition on Te- For the latter, setting u = u — S^g one has that 



- Aft + Vp = ^ eAwe.i - \/0^^, 

«e{/3.T,x} 

divu= ediv Wg^i + O 

ii - r = Ol 

P = J 



A5,g + £ 



in r2e, 

in 

on To, 
on Fjv- 



Thanks to the explicit form of Uq 



i 2:261 and to Proposition 15.31 one deduces that 



l|A5.g|lH-i(n,) = W^SMmn.^^ < fc {£i + (£ + £2)1 l0g£|^} , 
combining this with the results of Lemma l5.ll and using Theorem 14 . 5 1 one gets the desired result. 



5.1.8. Very weak estimates. We use here the framework of very weak solutions introduced 
above. The essential motivation comes from the lack of regularity of the averaged approximation 
(ucPg) across the interface Fq and the boundary layers' optimal cost in the x H^^ norm. The 
roughness Fe is contained inside the limiting domain 57i: we decompose our domain in three parts 
fl'i, Ce and fl2- 

Theorem 5.3. The full approximation {Ue,Ve) satisfies the error estimates: 



\u,-U, 



1UO2) + IIJ'^ ^'5llff-i(f2'iU£,U02) - 



where the constant k does not depend on e and | is a real strictly less than |. 

Proof. One gets thanks to Theorem 14. 7i very weak estimates on u := u^—U^ and p := p^—V^: 

l|u|lL2(o') + IIPllff-Mn') ^ + ll^ew,,, - Vei,,;!!^-!^^') + ll^div We,i||^2(n,) + 0(6^)} 

and 

I|u|1l2(£,) + \\p\\h-hc.) - ^^{Mhhc.) + I1p!1l2(£.)} • 

Thanks to Lemma |5. II and Theorem 15 .41 one has finally when gathering both inequalities: 
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< ke^' 



By a triangular inequality one obtains: 



(OiUi22) 
lff-i(^^iu£.U02) 



''^llL2(i2iun2) 



I_y-i(niua2) 



Next, setting vl := — = e(u| — (ui — g)) and p := 7'^'^ — V^ ^ eiPi — Pi), these variables 
solve: 



where again 



- Au + Vp = 
div u = 
u = eg(l - A,) 
U • T = eg(l - Ae) • t' 

ii = 0, 



in fliU Q2, 

on r2, 

on Tin U Tout,! U Tout, 2, 

on Ti UTo, 



dx2 



duo, 
dx2 



Using then again the very weak estimates of Theorem 14.81 one gets 

l|u|lL2(nin02) + llPilff-i(i2in02) ^^'e|lg(l - ^e)|lL2(r.„uro„t,iur2) + ^" ^\\drK\\H-^(T^) < ^'"'e'- 
We detail here only the second ter of the rhs above. If there exists h E Hq{Tn) s.t. 

d^h dX, 1 



(5.15) 



then ones has easily that 



As Fin U Tout,! are straight segments, an easy computation gives that if 







dh 


< 




dr 







h{x2) := ^l[o,e](a;2) + ^JfT^l[e,i](^2), Vx-2 e [0, 1] 



then h e i7g(rin U Tout,!) ^-^d it solves (|5.15p . An exphcit computation gives that 



dh 



L"(riv) 



which ends the proof. □ 

Remark 5.2. We are not allowed to apply the very weak framework to fii^^.' even for C°° 
obstacles, ^i,e does not satisfy uniformly wrt e the regularity property of definition Thus 
we applied the very weak estimates above the rough layer in this latter domain satisfying 
the regularity requirement of definition \4^.6\ uniformly in e. In the Ce zone we use the Poincare 
inequality to obtain the desired convergence rate. This explains why at last we obtain convergence 
results for the pressure terms in the H-'^{n[UCe^n2) norm which is smaller that the H ^(f2iUf22) 
norm used in the case of a flat sieve (cf. p. 50-52 in \12^ ). 

Here we consider the oscillating part of our approximation. We recall that Ue := Uq + eui and 
Pe '■ = Po + ^Pii and we set 

w^:=Ue-\ie, q^:^Ve-p^. 
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The functions Ve,gj, are explicit sums of all the correctors in (jS.lOp . In order to prove error 
estimates we need the following two results 

Proposition 5.5. // a periodic function p is harmonic on Z„oo,o o,i^d, on Z\_A^ao o,nd tends 
to zero when I2/2I goes to 00, then setting Pe =p(x/e), one has 

\\Pe\\H-^(n',) < WPAIh-HU;) ^ ^^^^ 

where the constant k is independent on e. 

Proof. We prove the result for ^I'l , the proof is the same for As p is periodic and harmonic 
in Zi^oo^ it is explicit in terms of Fourier series: 

It is thus decreasing exponentially fast towards 0. Then we solve the problem : find q s.t. 

-Aq^p in 

q = on {j/2 = 1}, (5.16) 

q is 1-periodic in the yi direction. 

Thanks to the exponential decrease of p it is easy to show that it belongs to D^'^(Zi.+oo)' and 
thus by the Lax-Milgram theorem, there exists a unique q € D^'2(Zi^+oo) solving (f5J6l) . One can 
even decompose q in Fourier modes and obtain again that it is an exponentially decreasing to zero 
at infinity. Then we set := q{x/e), and we have 

-Ax{e^q^)^p^ mQ,[. 

Given ip e HQ{n[), we aim at computing 

J((p) := / p^ipdx ^ I -Ax{e'^qe)ipdx = e"^ / \J ■ ^ x^-pdx, 
One has immediately because of the microscopic structure of q^ 



J{ip) / \/yq^-\/ipdx < eh\'^yq\\L2(z, 
Jn' 



the result follows writing that \\qe\\H-Hn[) = ^^PveHl,{n[}iJiv)/\\v\\HUn[)) ^ 
For the vertical correctors of pressure terms one has in the same way: 

Proposition 5.6. For a given G VF°'^(n) s.t. a €]0, 1[, setting 0^ (^{x/e) one has that 

where the constant k is independent on e. 

Proof. We restrict ourselves to the case of ft'i again. We solve at the microsopic level: 

i-Aq = e, inM+x]l,+oo[, 
\q = on {0}x]l,+oo[UlR+ X {!}. 

With arguments similar to those of the proof of Proposition lA.ll one can show that if 9 is in 
W^~"^'^(M+x]l, +oo[) with S g] — 1; 1[ then there exists a unique solution q € ^^^'^(R+x]!, +oo[) 
solving (|5.17p . An easy computation shows that if 6* e 11^°'^ (11) then 9 e W^^'i{Il), which implies 
setting S ^ a — 1 the existence of a solution q G M^2-i(n) provided that a €]0,2[. As, by the 
definition of 9, a e]0, 1[, we restrict ourselves to solutions q G W^!^;^(n) with a s]0, 1[. Again we 
set ge = q{x/e) which means that 

-A,{e\,) = 0, in Q[. 
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Given a test function Lp £ Hq{D.'i), we aim at computing 



Because of the microscopic structure of one has again 

J{ip) = <^ "^yQe ■ Vx{i^ip)dx < e|lV'llw/i.=o(o;)l|Vj;(7e|li2(o-)l!^|lHi(n;)' 
passing from the macro to the micro scale we have 



Ji 



V,qe|li.(f,,) < ( / / |Vq|V"-'rf2/ sup 



pG-B(0, 



by similar arguments as in Lemma 15.11 Again the result follows writing that ||9e|li/-i(n' ) 
Theorem 5.4. The rapidly oscillating rest (Uc — Ug,Pe ~Pt) satisfies 

where the constant k is independent on e. 
Proof. Because is explicit and reads : 



duo. 



dxo 



a direct computation of the norm gives that 



1v,|1l2(o^) <ek 



/3e-/3 



L2(Oj) 



3 



L2(0,) 



- e >ff — x\ 



<=I|w^,«IIl2(o,) + I1W||l2(o^) < fee 



We use again the decomposition of in subdomains f2i,£e and 0,2- The t/i-periodic pressures 
TTe, zue,r]f — 77 and /ie fulfill hypotheses of Proposition l5.5l the vertical correctors 9i for i € {/3, T, x} 
satisfy hypotheses of Proposition 15.61 one then concludes 



where the pressure correctors as 5 and e terms are implicitly treated by a direct estimates of the 
norm. In we use the dual estimate (|4.5p based on the Poincare inequality, to get 



Combining Theorems 15.31 and 15.41 above, one gets the proof of Theorem 13. II 



5.1.9. Implicit interface conditions. We start with the horizontal velocity. We call u^^ 
(resp d2U^i) the values above and below Fq. The first order interface condition derived above on 
Fq reads: 



u 



I dx2 



dx2 
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=± 



ei 



bo] 



assembling together normal derivatives of the velocity on both sides and because 0x2^^ i = 0, one 
has also : 



which finally gives 



dx2 
dx2 



ei 



ei 



-S-X2e2, 

m 

[Po] = 

-1^X2^2, 

m 



and 



uf ■ ei 



Setting := Uq + eui and because Uq = on Fq, one has also 



dx2 



+ O(e^), and 



One recovers a slip velocity condition in the main artery and a new discontinuous relationship 
between the horizontal components of the velocity at the interface. 

For the vertical velocity, thanks to the continuity of ^2 across Fq, one has that 

bo] _ ([cruo,po] ■ 62,62) _ ([cru,,pj • 62,62) 



n,2 



n,2 



Ul.2 



■0(e), 



this in turn gives the implicit interface condition : 



■ n = -^([ctu,,pJ •n,n) + 0(e^). 

m 

5.2. The case of an aneurysmal sac . When e goes to 0, the limit solution (uo,po) is 
explicit (we set Pout.i = in p.6p ): 



Vio{x) = -^{l - X2)x2eiln^, ^x e fl 

Po{x) = pin(l - a;i)loi +Po lf22; 

where is any real constant. Following the same lines as in Theorem 15 . II one obtains 

Theorem 5.5. For every fixed e, there exists a unique solution {u^,pe) G x L^(f2g) of 

the problem (j3.ip . Moreover, one has 



l|Uf - uo|lHi(n,)2 + be -PollL2(jy) + V^\W- Pa\\L--{c,) < fcVe 

where the constant k depends on Pq but not on e. 

5.2.1. First order approximation. Due to the presence of three kind of errors above, we 
construct a full boundary layer approximation {U^,Ve) exactly as in (j5.10p . One has to make few 
minor changes in the definition of {WejZf) that are left to the reader. The only difference stands 
in the pressure jump: 

[Po] =pj(a:^i,0) -Po , 

where p^ is the constant pressure not yet fixed. The first order macroscopic corrector (ui,pi) 
should satisfy 



-Aui + Vpi = 
div Ui = 
ui = 
ui • r = 0, 

Pi =0, 



Ul 



dx' 



-/3 



=± 
T 
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[Po] = 



in ill U f22, 
in r^i U ^2, 

on Fi U F2 U Fout,2, 

on Fin U Fout.i, 
on Fo^. 



(5.18) 



As we impose the velocity on every edge of ^1,2 there is a compatibihty condition between the 
Dirichlet data and the divergence free condition reading 

/ divui(ia;= / Ui-nda— ui-n(ia;i=0, 
Jn2 JdQ.2 Jto 

and this precisely identifies the pressure giving 

|rob,7 = / P^(xi,o)dxi. (5.19) 

The first order constants arc fixed in the definition of (U^jVe)- Even if Pq is now well defined, 
the first and second order pressures pi and p2 are again computed in fl2 up to a constant. This 
is why we still need norms on a quotient space L^(J72)/K for the pressure in $72- Following the 
same lines as in the section above but taking into account the pressures in up to a constant as 
in the proof of Theorem 15.51 one proves Theorem 13.21 

6. Numerical validation. We present in this section a numerical validation in the case of a 
collateral artery, as one obtains similar results in the case of an aneurysm we do not display these 
results. We solve numerically problem (|3.1|) in 2D, for various values of e. For each e, we confront 
the corresponding numerical quantities with the information provided by the homogenized first- 
order explicit approximation : velocity profiles, pressure, fiow-ratc. Numerical errors estimates 
are computed with respect to the different norms evaluated above in a theoretical manner. 

We do not include in these sections approximations based on the implicit interface conditions 
presented in p.4|) : this will be done in a forthcoming work that investigates new theoretical and 
numerical questions that these conditions pose. 

6.1. Discretizing the rough solution {ug,p^). The domain ft^ is discretized for e e]0, 1] 
using a triangulation. To discretize the velocity-pressure variables, a (P2, Pi) finite element basis is 
chosen. Because of the presence of microscopic perturbations, when solving the Stokes equations, 
the penalty method gave instabilities. For this reason we opted for the Uzawa conjugate gradient 
solver (see p. 178 in [T7], and references there in). The code is written in the f reef em++ languag^. 
On the boundary we impose the following data : p-m = 2,pout,i = 0,pout,2 = —1- In order to 




Figure 6.1. Direct computation u^^i (left) Uf^2 (middle) and (right) 

improve accuracy of the direct simulations we use mesh adaptation iterations as described p. 96- 
97 in |17| : using the hessian matrices of components of Ue , one defines a metric that modifies the 
mesh (see fig. 16.71 (middle) for a final shape of the mesh). 

6.2. The microscopic cell problems. Using the same numerical tools, we solve the micro- 
scopic problems (j5.2p . (|5.3p and (|5.4p . These are defined on the infinite perforated strip Z: one is 
forced to truncate the domain and works on Z-l^l with L > large. We impose boundary data 
at the top and the bottom of Z^^ ^ namely 

^2(2/) = T2(j/) = 0, X2(y) = -1 on{2/e]0,l[xRs.t. y2 = ±i}. 



^http : //mm . f reef em ■ org/ f f 
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and we let natural boundary conditions on the other components. When L goes to infinity it is 

proved in |23| that the solutions of the truncated problem defined on Z^l l converge exponentially 

=± 

with respect to L to the solution of the unbounded problem. We compute numerical values of Pi 

=± _ 
and T]^ and the pressure drop [77]. If Js is a sphere of radius 3/16 in a period of size 1 centered 

at (1/2, 1/4) the numerical computations provide values listed in table WA\ One can notice that 



constants 


values 


constants 


values 


Pi 


-0.377928 




-0.122114 




-0.000371269 


% 


0.121744 


m 


27.9435 







Table 6.1 
Homogenized numerical constants 



contrary to the resistive matrix of [5] the tangential part of the coefficient are negative. This is 
due to the fact that the obstacles lie above the interface in the main flow. The horizontal first 
order slip velocity is thus negative (see below). 




Figure 6.2. Velocity vectors for micrscopic problems: (3 (left), T (middle) and x ('"S^iJ 




Figure 6.3. Pressure for micrscopic problems: n (left), ii(middle) and rj (right) 



6.3. Explicit first order problem. Wc solve problem p.3p on triangulations of 51i and fi2- 
Because of the discontinuity of the Dirichlet data at the corners O and x, the solution (ui,pi) 
does not belong to H^(r2i U ^,2) x L^{ili U fl2)- Indeed a pressure singularity occurs at O and x: 
refining the triangulation at the corners one get a point-wise explosion of the pressure near O and 
X. We add then the zcroth order explicit poiseuille profile to obtain a numerical approximation 
of (ue,Pg). For e = 0.25, we display in fig. 16.41 velocitv components and pressure projected on 
in order to be compared to (u,;,^,;) in the next paragraphs. Since the pressure is not bounded 

29 



Figure 6.4. Explicit first order approximation u^^i (left), itE_2 (middle) and (right) 



(the numerical value is very high in a very small neighborhood of O and x we display the "regular 
part" : we cut-off the pressure function near the corners for visualisation purposes only. 

6.4. Comparisons and error estimates. In fig. 16.51 we display the horizontal (left) ve- 
locity profile above and below the obstacles for the direct solution Ue^i and our approximation 
Ue,i- In the middle we show the normal velocity in the same framework. On the right of the same 
figure, we plot various values of e on the x-axis and the flow-rate through Tq on the y-axis. One 
observes that the asymptotic expansion gives the first order approximation of the flow-rate with 
respect to e near e = which was expected. One notices also that the actual rough flow-rate 
behaves as a square- root of e. This seems difficult to prove using averaged interface conditions 
only [51 [7]. In fig. 16. 6[ we plot numerical error estimates for the zero order approximation (uo,po) 
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Figure 6.5. Velocity profiles for e = 0.125 close to Tq: horizontal (left) and vertical (middle). Rough versus 
homogenized flow-rates (right) wrt e 

and for our explicit first order averaged approximation (ue,p^) wrt to the direct solution (uc,pe). 
Left we display the L^(f2c) error for velocity vectors. On the right, we compute the pressure error 
estimates in the H^^{il[ U £c U fl2) norm: for po (rcsp. p^) we solve numerically 

- Aq pe - Po (resp. pj, in n[ U £e U ri2, 
q = on af^; UPU9r22, 

for each e then we display ||V(7||^2fQ/ u£ un V recovers theoretical claims of Theorems 15 . II and 

In fig. 16.71 middle and right we display the meshes used for a single value of e = 0.25 and for 
the computations of (ui,pi). On the left we display the mesh size h used for the direct simulations 
with respect to e. 
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1 0.01 



Figure 6.6. Numerical error estimates : velocity (left) and pressure (right) wrt e 





Figure 6.7. The mesh size h for direct computations (left), a direct mesh for e = 0.25 (middle) and the mesh 
used to compute (ui,pi) (right) 




Appendix A. Well-posedness in weighted Sobolev spaces of the vertical correctors. 

Given the data {i,h), we study the problem: find {w,6) e W;!,'^(n) x W^'^(n) solving 



(A.l) 



Remark: 2. For the specific type of mixed boundary conditions set on dll, it is not possible 
to use neither the Fourier transform in the vertical direction nor the Laplace transform in the 
horizontal direction in order to derive results established below. To our knowledge there are few 
results in the literature for unbounded domains with mixed boundary conditions on noncompact 
boundaries. 

Theorem A.l. // the real a is such that |a| < 1 and if 

i e'Wl'^{D\JN), h€W~^''^{N), 



there exists a unique solution (w, 0) e W;!j,'^(n) x py^'^(n) solving problem (|A.ip . 

Before giving the proof of the theorem, we need to prove two intermediate propositions. For 
this sake we define 

X„ := {v e W;!,^2(n) s.t. V = on D, v • r = on iV}, := W^"^ {li) . 

First we show that the divergence operator is surjective from into Yq 

Proposition 7. For any given function q ^ there exists a vector function v G such 

that 

divv = g, and \^\xa < ^(n, a)||(j||y^ 
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where the constant k depends only on the geometry of the domain, and on a. 
Proof. We define a sequence of annular domains covering 11 

C„ := {y e n s.t. if y = (r, ^) tlicn re]2"-\2"[}, n>l, Cq ~ B{0,1) DU. 

We decompose q as q ^ X^^^o?"' '^ith qn ■= 1^c„- On each. C„ we solve the problem: find 
G JCa^7i s.t. div V71 = qn and Iv^^jj^^ 

< k{Cn , Q!)|'Zn|y^ , where 

X„,„ := {v e X„ s.t. V = on {lyl = 2"-i} U {\y\ = 2"} U (C„ n D), 
and V • T = on (C„ n A^)} . 

But to solve the latter equation in a weak sense means 

/ div Vn ■ OJ rdrd0 = / qn-ojrdrd9, Vw G W^'^(C„), 

JCn JC„ 

making the change of variables: (f = r/2"^^,6) and setting 

v„(f , 9) :== v(2"-i7^ d), Ur, 0) := q{2"-^f, 9), tZ.(f , 9) := L.(2"-if, 9), 
the problem becomes: find v G Xo.i defined on Ci s.t. 



divv„-(I)f(ird6l = 2""M q-Cjrdrd9, Vw G i^(Ci), 

the test space is defined on a compact fixed domain Ci, weighted Sobolev spaces coincide with 
the classical ones as soon as the weight is strictly positive and bounded. In this framework the 
operator div : Xq i — > Yq i is surjective thanks to Lemma 4.9 p. 181 in |13j . Thus there exists 
v„ G Xo.i s.t. divv„ — q„2"~^. Note that there is no need of a compatibility condition on the 
integral of g„ as in Lemma 3.1 chap. Ill in [TS] because v„ • n 7^ on n Ci. Moreover one has 
that 

I^kIh^Ci) - ''^(^l' 0)||2" "''?n||i2(c'i)' 

where k depends only on the geometry of Ci and is thus independent on n. Turning back to the 
original variables (r, 9) one has then that div v„ = g„ and 



/ \S/Vn\^rdrd9 <k{Ci,0) [ \qn\^rdrd9. 
Jc„ JC„ 



In order to recover the global weighted norm of q in W^^'^ljl), we multiply the inequality by 
22Q(n+i) both sides; we use that for r G [2"^^, 2"], p := (1 + r^)2 can be estimated as 
22a(n-i) < < 22"("+i) giving finally 



<fc(Ci,0)22"("+i) / qldy<2'^k{M) qlp^'dy. 

One defines v := v„lc„, because of the boundary conditions imposed on each of the C„, v is 
continuous on 11 and thus belongs to ~W]^^{J1). This gives the result. □ 
We lift problem (|A.ip by subtracting to w a function 7?,(w) satisfying: 

7e(w) G 7^(w) = f on £> U A^. 

Such a lift exists (cf p. 249 [16] for an explicit form of 7?.(w)). We correct the divergence of 72.(w) 
by setting : 

S{vf) G Xa S.t. div5(v^f) -div(7^(v^r)) and i|5(w) ||^i^,2(jj^ < fc||7^(w)||^l^,2(JJ^, 
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which is possible thanks to Proposition [7l The new variables (w := w — 7?.(w) — S{w),9) solve 
the homogeneous problem: 



- Aw + V6' = A7^(w) + A5(w) in n, 

divw = in H, 

w = on D, 

w ■ T — Q, and 9 = h on N. 



(A.2) 



No we claim that (F, G) 

wl'\u) X <'^(n) s.t. 



- (p^-w, p°'6) solve in an equivalent way the problem: find (F, G) in 
^„ F + BlG = p"(A7^(w) + A5(w)) 

m ii, 

(A.3) 



where 



F = 
F • T = 

G = 

1 



in n, 
on D 
on N 



1 



AF-2p"VF-V p"A— F, and S^F := divF + p"V — • F. (A.4) 



1 



Indeed, p € C°°(n) thus if {w,9) solves (jA.2p in the distributionnal sense, then equivalently by 
its definition the pair (F, G) solves (|A.3[) also in the distributionnal sense. Uniqueness is insured 
thanks to the onto mapping between W^^^^n) x W°''^(n) and Wj'^(n) x W^^'^in) (cf Theorem 
1.3 p. 243 in [l^): if {w,9) is a unique solution of (|A.2[) then so is (F, G) for system (|A.3p 
and vice versa. The boundary conditions match between both problems by similar onto trace 

mappings. Note that the rhs in (jA.Sp belongs to 'Wq^'^{II) and the boundary data to Wq ^'^(iV). 
We associate to (jA.3[) the corresponding variational setting, namely we define: 

• the velocity /pressure test space is Xq x Yq, 

• the bi-continuous (resp continuous) forms aa,ba (rcsp. la) read 



a«(F,V) = (VF, VV) - 2 p"VFV 



p"A 



1 



F,V 



VF,Ve Xo, 



6a(F,( 



divl 



F,Q , VFeXcVQeFo, 



[ /„(V) = (p"(A7^(w) + A5(w)),V)n - (p"/i,V- n)^, Wg Xo. 

• the variational problem: the problem (jA.3p can then be restated in an equivalent way: 
find (F,G) e XqxYo s.t. 

ra„(F,V) +fe„(V,G) = /„(V) WeXo, 
|&„(F,(Q) =0 VQeFo- 

We denote by Aa ■ Xq — > X'q the operator s.t. 

a„(F,V) =<^F,V>x;,,x„, VV e Xq. 

The well-posedness of problem (|A.3[) is equivalent to two conditions (Theorem A. 56 p. 474 |13)): 
(i) PAa ■ ker(;BQ) kei^Ba)' is an isomorphism 
(a) Ba ■ Xq — >■ Yq IS surjective 
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where P is the restriction of Aa to the kernel of Ba- Here we prove that these conditions are 
actually fulfilled. 

Proposition A.l. If\a\ < 1 thenAa satisfies condition (i), whereas B a satisfies {ii) without 
any restrictions on a 

Proof. We prove at first that condition (i) is satisfied by showing that is a coercive bi- 
linear form. For every vector F in Xq, one has after integration by parts of the second term in 
the definition of Or,: 



a„(F,F) = I ai^-:^F,F 



IVFI 



an 



a -div ^ + ^ - 



(^ + 2)|yp 



note that the boundary term on the first line above vanishes on N because (y • n) = 0, though 
Fi 7^ on this part of the boundary. Note also that this integration by part is justified for functions 
in Xq. We use optimal Poincarc-Wirtinger estimates already presented in the proof of Theorem 
5.3 p. 20 in [27]: 







2 


F 


2 




dy < 




< 








P 


L2(n) 



< |VF| 



L2(n)- 



Note that these Poincare-Wirtinger estimates are possible because of the homogeneous Dirichlet 
conditions on D: they give stronger weights than the corresponding logarithmic weighted Hardy 
estimates available in the whole [3]. Finally one has 



(^F,F)n > (l-a2)|VF|^, 



which implies coercivity of the operator if |q;| < 1. Note that this result (also valid in the scalar 
case) improves Lemma 4.3 in [S]. This is essentially due to the integration by parts performed on 
the term (VFy/p^, F)n which avoids estimating this term separately from the others. 
We focus on the condition (ii). For all g € we look for F S Xq s.t. 



Bo^W^q, and ||F||^^ </c||<j||^o,2 



(n)' 



but this is equivalent to solve 



div — = ^ 



li q E Yq then q/p" E W2'^(n) and by Proposition [7] there exists v E X^ such that 



w°'^(n) 

1,2/ 



Set F := p"v thanks to the isomorphism between W;|,'^(n) and Wq' (H) there exists a constant 
s.t. 



|F|w-(n) 



< ^i|v|wj,'2(n) ^ fcifc(Ci,a) 



^•'l9l<'"(n)- 



□ 

Proof. [Proof of Theorem lA.lj Thanks to the equivalence between well-poscdness and con- 
ditions (i) and {ii) one concludes the existence and uniqueness of a pair (F, G) solving problem 
(|A.3p . Moreover one has the a priori estimates : 



Iwi-^(n) 



IMwrin) < k' (^|lp"A(7e(w) +5(w))||^-i,.(n) + \\p"h\\ 



{N) 
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they are obtained similarly to those of Theorem 2.34 p. 100 in [13]. The isomorphism between 
weighted spaces mentioned above and the equivalence of problems (|A.3p and (|A.2|) gives existence 
and uniqueness of (w,0) solving problem (jA.2[) and a priori estimates 



l|w|lwi-^n) + PWw^^n) < fc" \ !|A(7e(w) +5(w))||^-i..(n) + 



< fc'" |||7e(w) +5(w)||^i,2(n) + \\h\\^_.^,^^^ 

This gives existence and uniqueness of (w, 6) and due to the continuity of the lifts 72.(w) and 5(w) 
with respect to the data, one easily proves that 



l|w!lw-(n) + ll^llM.r(n) < |l|f|lw|-(.uA.) + ""^t-i^w 
which ends the proof □ 

Proof. [Proof of Theorem I5.2j We use Theorem lA.ll to prove results for {wi,di) for i e 
{/3, Y, %, The data for these problems tends exponentially fast to zero: in terms of weights, 

- 2 -i 2 

the Dirichlet (resp Neumann) data is thus compatible with Wa {D) (resp. Wq ^' (N)) for any 
real a. Setting f = (3X — f3 (resp. f = TA — T and f = xA — x ) hi (jA.ip is equivalent to the first 
(resp. second and third) problem in (|5.5p . Applying Theorem I A . 1 1 gives then the claims. □ 

Appendix B. Periodic boundary layers: proofs of Propositions [2] and [3] and of 
Corollary 15.11 

Proof, [of Proposition [2] We start by lifting the non-homogeneous Dirichlet boundary con- 
dition : we set TZ{(3) :— ?/2<?!'(y2)eilz+ and ^ := (3 + 7?.(/3), this is still a divergence free vector. 
Now, there exists a unique solution s.t. 



V/3 

Indeed, in the space of divergence free D^''^{Z) functions vanishing on P, the gradient norm is 
a norm (via Wirtinger estimates), one proves existence and uniqueness of /3 in D^''^{Z) by the 
Lax-Milgram Theorem. 

We apply Lemma 3.4 and Proposition 3.5 of [TSl in order to recover the Lf^^{Z) pressure 
solving: 

-A^ + Vt: = -A7^(/3), 

and this gives existence and uniqueness of 0,7r) e DI'^{Z) X ifoc(Z). On the interface located 
above (resp. below) the obstacle J7s we apply the Fourier decomposition in modes as in Theorem 
3 p. 10 [24]. One obtains the exponential convergence towards the zero modes of (/3, tt) in an 
explicit way. To derive relationships between constant values at infinity, one has 

(i) by the divergence free condition that f^2{'^) = 1^2 (t) = for all v > 1/2. p and 7 < 0. 

(ii) integrating the first equation of (|5.2p in every transverse section {1/2 = S} which does not 
cross the obstacle Js gives 

/ /■ o , \ f ^'/^i 9tt ^ 

dyi \J{y,=s} J J{y2=5} m 9yi 

by j/i-periodicity. This implies that /3{d) is an affine function of d. As the gradient rapidly 
goes to zero, the linear part is zero, we conclude that only the constant remains : thus 
Pii^) = ^i(+oo) for S > y2,p, and '^^{1^) = '^^^{0) for v < 0. 

(iii) Set G := y2'^z+^i and F := 0, they satisfy: 

= — (5y^ei in Z, 

. ; (B.l) 

m Z, 




we test the first equation in (|5.2p by G and the first equation in (jB.ip by /3, then we 
integrate on : 

(AG - VF, /3) - (A/3 - Vtt, G) = (<5sei, /3) = /3i(2/i, 0)dyi = ^i(O) 

Jo 

= (ctcf • n, ^)9z,,^ - (cr/3,7r • n, G)az„^-, (f^cF • n, /3)p + /S^ (;^) - {ap,^ ■ n, G)p 
= -(9n (y2ei),y2ei)p +^1(1^) + (a/g.^r • n,/3)p 
= -(9n (y2ei),y2ei)p +^1(1^) + |V/3|i2(2), 

where we neglected exponentially small terms on {j/2 = i^} and {2/2 = 7}- Now we explicit 
the physical meaning of the constant Q := ((?„ (2/261), 2/261 )p 

S+ (<9„ (2/261), 2/2ei){y,=^}u{y2=7} = (^(2/261), J/2ei)z„^ + |V(2/26i)||2(2_^^)4, 

which in turn gives : 

where \Js\ is the volume of the obstacle Js- The quantity Q represents the volume of fluid 
missing due to the presence of the obstacle Js above the limit interface S. If we were to 
consider a straight channel without a collateral artery but a roughness below the fictitious 
interface, Q would be a positive number equal to the volume of fluid present below E. 
Computations above are formal and can be rigorously derived by regularizing the obstacle 
Js and then working on regular functions in order to obtain results stated above. None 
of the final quantities depending on second order derivatives, passing to the limit wrt to 
the regularization parameter, extends results above to Lipschitz obstacles. 

□ 

Proof, [of Proposition |3] The existence and uniqueness part follows exactly the same lines 
as in Proposition [2l the exponential convergence is also proved the same way. We detail only 
relationships between horizontal averages. 

• For T2 one uses as for the divergence free condition together with the boundary 
condition at infinity to obtain that T2(2/2) = for all 2/2 in 1R\]0, 2/2, p[- 

• Testing the first equation in (j5.3p by T and integrating on Z^^^ one gets when passing to 
the limit wrt v — > 00, 7 — > —00 that: 

Ti(0) = ||VT||^2(;,). 

Testing the same equation again but with G = 2/261 and integrating on gives: 
Ti(zy) = Ti(7) + (cty.ct • n, 2/2ei)p, Vi^ > 2/2,p, V7 < 0. 

We compute: 

(div(Tr,tu,/3)z„,^ - (div (7/3,^, T)z,^ ^, = -/3i(0) 
=(o"r,ro ■ n, -2/2ei)p + (crx.ro ■ n, /3){,y2=iy}u{y2=7} ~ {'^P.t ' If ){t^2=i'}u{y2=7} 
which gives after passing to the limit wrt v and 7 

/3i(0) ^ (cTT.ra • n, j/26i)p. 

• Testing the first equation in (|5.3p against 2/2 Iz- and integrating on Z^^^ one obtains easily 
that 

Ti(2/2) = Ti(0), V2/2<0. 

Putting together equalities obtained above one concludes the proof. 
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n 

Proof, [of Corollary 15. 1[ Setting again := f3 + 2/261 lz+ and writing 



(Ax - Vr/, ^)z,^„ + (A/3 - Vtt, x)z.„. = 

(f^x.r) ■ n'/3){a2='5}u{y2='4 ^ (^/a.TT ■ n, x){a2=5}u{a2=>'}- 
When passing to the limits S 00 and — > — cx3 in the last expression, one obtains that: 

now because — >■ and X2 ^ ^1; one gets the desired result at infinity. As the pressure tt is 
harmonic in Z, the average 7f((5) is zero in R_U]?/2,P7 +oo[. The same proof holds for zu. □ 
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